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Abstract 

The degree by which a function can be differentiated need not be re- 
stricted to integer values. Usually most of the field equations of physics 
arc taken to be second order, curiosity asks what happens if this is only 
approximately the case and the field equations are nearly second order. 
For Robert son- Walker cosmology there is a simple fractional modifica- 
tion of the Friedman and conservation equations. In general fractional 
gravitational equations similar to Einstein's are hard to define as this 
requires fractional derivative geometry. What fractional derivative ge- 
ometry might entail is briefly looked at and it turns out that even 
asking very simple questions in two dimensions leads to ambiguous or 
intractable results. A two dimensional line element which depends on 
the Gamma- function is looked at. 
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1 Introduction 



1.1 Motivation 

The field equations of fundamental physics are usually taken to be second 
order. Like any other property of a physical theory this should be sub- 
jected to experimental and observational tests to see what the experimental 
bounds are. A method of examining this is by using fractional derivatives to 
investigate the properties of differential equations which are almost second 
order. Electromagnetism might produce the best tests of how near to second 
order governing field equations need to be, but here an attempt is made to 
see what modification of gravity theory can be constructed using fractional 
derivatives. 

1.2 History 

In 1695 L'Hopital wrote to Leibniz asking him what happened if the num- 
ber of times a function was differentiated was not an integer but rather ^, 
from this early beginning the subject of fractional derivatives was born. Al- 
though the derivatives are called fractional derivatives they can take any real 
or sometimes complex value. Recent textbooks include [-1, 7]. Recent ap- 
plications to physics include construction of a fraction Schrodinger equation 
[6, 4] and some properties of field theories [1, 11]. 

1.3 Methodology 

There are two distinct methods of approaching what fractional derivative 
cosmology could be. The simplest is last step modification in which Ein- 
stein's field equations for a given geometric configuration are replaced with 
analogous fractional field equations, in other words da — > after the field 
equations for a specific geometry have been derived. The fundamentalist 
methodology is first step modification in which one starts by constructing 
fractional derivative geometry. The problem with the last step approach is 
that it probably only gives consistent answers for geometric configurations 
that are expressed in rectilinear coordinates; the problem with the first 
step approach is that the whole of geometry has to be rethought through, 
even things like linear coordinate transformations have to be replaced by 
fractional ones, perhaps quadratic forms by fractional forms and so on; in- 
termediate step appraoches seem to have the disadvantages of both of the 
above approaches. When using fractional derivatives it is not always the 
case that differentiating a constant gives zero see Podlubny [7]p.80, so it 
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is not clear what kind of fractional derivative to use, here a pragmatic ap- 
proach is adopted: if a fractional derivative can be made to work then it is 
used. 

1.4 Outline &z Conventions 

In §2 an attempt is made to see how the newtonian ^ potential is modified 
when the d'Almbertian is replaced at the 'last step' by fractional deriva- 
tives, it is found that the 'last step' assumption is not necessarily correct as 
spherically flat spacetime need not take its familiar form. In §3 Robertson- 
Walker cosmology is looked at in the last step approach. In §4 briefly looks 
at what fractional derivative geometry might entail. In §5 is the conclusion, 
p is used for degree of fractional differentiation, the familiar non-fractional 
differential equations are recovered when p = 1, V,n,Ua are the pressure, 
density and co-moving vector of a fluid. Calculations were carried out using 
maple9 / grtensorll [5] . 

2 Newtonian Gravity 

Vacuum Newtonian gravity is taken to be governed by the D'Almbertian 
operator acting on a scalar function 

□</> = 0, (1) 

Applying the last step fractional derivative substitution approach for spher- 
ical symmetry gives 

p2p^ + ^P^V = (2) 
r 

numerically solving this for values close to one produces potentials as illus- 
trated in figure one. For p close to one the potential can be made arbitarly 
close to the non-fractional case. Rates of decay of potentials have been dis- 
cussed at lenght in [8]. At first sight one might imagine that all that now 
needs doing is using either expansions or numerical analysis to compare this 
result with observations: but there is there is a serious problem with equa- 
tion 2, the ^ term assumes a standard spherically symmetric flat background 
spacetime, however for fractional derivatives this assumption does not nec- 
essarily hold, the two-sphere in spherical coordinates is not necessarily the 
familiar one as this requires assumptions of standard geometry which might 
not hold in fractional geometry, see §4.2. 
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3 Fractional Robertson- Walker Cosmology 

3.1 Friedman Equation 

For the moment we retain the standard quadratic form of the Robertson- 
Walker hne element 

where fc = 0, ±1. The stress is taken to be that of a perfect fluid 

Tab = {fJ^ + r)UaUb+Vgab, Ua = (1,0,0,0) , (4) 

which for present purposes contains no derivatives, this might change if 
the Clebsch representation of the comoving vector field Ua is used. For 
general relative the dynamics are governed by the Friedman equation and 
the first conservation equation UaT"'!^ = replacing partial derivatives with 
fractional derivatives at the last moment the Friedman and conservation 
equations become 

3 [fc + (Pfa)^j = Kfia^, a^V^V = [a^ + V)] , (5) 

respectively. For the simpest non-fraction p = 1 case take both k,V = 
then the conservation equation integrates to 

^ = -^, (6) 
and then the Friedman equation integrates to give 

For the =p 1 cases powers of t are replaced be trigonometric or hyperbolic 
functions. 

3.2 Naive Approach 

A naive generalization of the above non- fractional approach is just to re- 
place the time derivative in the Friedman equation with a fractional time 
derivative. The Riemann-Liouville fractional derivative with lower terminal 
at 0, after correcting a typo [ ] page 310 is 
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Absorbing constants into C gives the scale factor and density 

a = C7tf, /i = C-V2P, (9) 

which is what might have been anticipated. Tlicre is a problem with the 
above in integrating the conservation equation 6 it was assumed that differ- 
entiating a constant gives zero but 8 gives 

which is non-zero for non-zero B. 

3.3 Compensating Pressure. 

To rectify this a compensating pressure can be used, substituting 9 back in 
to 5 a pressure which satisfies the conservation equation is 

AC^ is plotted in figure two. This shows that the compensating pressure is 
smaller than the density. 

3.4 7-equation of state. 

The solution 11 is an example of a 7-equation of state for which 

p=(7-l)/i, (12) 
for 117 = AC^ -|- 1. To investigate if 11 can be generalized choose 12 and 

a = At"", n = Ce, (13) 
then the Friedman equation restricts the values of n and C 

„=_2p, C=5^I(!^, (14) 

KT[m+ 1 — p}^ 



and the conservation equation places the restriction 

(7 - 1) r(l - 2p) _ r(3m -2p+l) 
7 r(l - 3p) ~ r(3m-3p+l)' 



(15) 



because m occurs inside a F-function it is hard to access the use of this 
solution. Generalizing 13 with V = the conservation equation gives 
n = r so that the results are the same. 
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4 Fractional Derivative Geometry 



4.1 Fractional alterations 

A generalized Christoffel symbol can be taken to be 

nc = \9"^ {^cQM + V^Qcd - V'gbc) , (16) 

where Dc is an operator with one spacetime index, in particular it can be 
a vector field, or the familiar partial differential operator, or a fractional 
derivative. A generalized Riemann tensor can be taken to be 

GRbcd = "^cTdb - T^d^ch + ^tf^db ~ ^2/^4, (17) 

where for simplicity the operator Vc is chosen to be the same one as in the 
generalized christoffel symbol. The fractional quadratic form could be 

ds^P = [dx'Y {dx^Y , (18) 

for p = the metric requires 2 x a indices, otherwise the number of 

indices is not integer, the meaning of a metric with non-integer number of 
indices is not clear. Forms of the form 18 can occur in Finsler geometry [2], 
however in Finsler geometry there are still linear coordinate transformations 
rather than fractional ones of the form 24. 



4.2 Fractional Flat Space 

It is now possible to ask what is fractional rectilinear flat spacetime, for 
simplicity working in two dimensions with the usual quadratic form 

ds^ = -N{t, xfdt^ + F{t, xfdx'^, (19) 

if the generalized Christoffel symbol 16 vanishes then the generalized Rie- 
mann tensor 17 will also vanish, looking at one component of the generalized 
Christoffel symbol 

^tt = ^^tN' (20) 

this will not vanish unless either = in which case 19 degenerates or 
caputso derivatives are used which results in the familar two dimensional 
rectilinear flat space with N = F = 1 . The next question is what is fractional 
spherical flat spacetime. In two dimensions this question is what is the 
analog of the following: start with line element 

ds^ = dx^ + dy^, (21) 
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and then perform the transformation 



x = rsm{0), y = rcos{6), (22) 

to give the hne element 

ds'^ = dr^ + r^d9^. (23) 
For the fractional case start with 21 and use coordinate transformations 

dx' = Vfx'dx'\ (24) 

choosing 22 to be replaced by 

x = r^sin(0), y = r^cos(0), (25) 

gives the line element 

ds^ = T{p + ifdr^ + 2T{p + l)rP cos{^)drde + r'^'PdO'^. (26) 

Altering the power of the trignometric parts of 25 does not seem to remove 
the cross term but rather makes things more complicated. Maple/grtensorll 
[5] works out the standard, non-fractional derivative curvature to be 

Rab = ^Rgab, Gab = 0, RiemSq = R^, RicciSq = ^R^. (27) 

Like all two dimensional line elements this solution obeys the vacuum Ein- 
stein equations. For fixed constant p it has constant curvature. Using the 
method of signature constants [9] the signature of 26 gives the same curva- 
ture in any combination. Two dimensional line elements have been studies 
by Witten [10]. Whether the line element 26 is explicitly flat using the 
curvature 17 leads to a calculation which has so far proved intractable. 

4.3 Fractional geodesies. 

Instead of starting with 16 and 17 one could attempt a more systematic 
approach by starting with a point particle lagrangian 



and then varying the fractional alteration to form a geodesic equation. 28 
can also be thought of as the integral form of an element of arc, compare 
equation one [ ']. Once one has a geodesic equation one has a connection. 
From the connection the commutation of covariant derivatives will give cur- 
vature. The problem with this is that it is not clear what the right hand side 
of 28 means, in principle it is not necessary to know explicity what 28 means 
only to know what the variation of 28 is, this is unclear because terms such 
as (^)^ occur and there appears to be no way of defining them. The same 
problems seem to occur starting from two other forms of the point particle 
lagrangian such as the second-order form and the Hamiltonian form. 

5 Conclusion 

It is possible to produce fractional generalizations of Newtonian mechanics 
and Friedman-Robertson- Walker cosmology by replacing partial derivatives 
by fractional derivatives in familiar equations. In cosmology the degree of 
fraction differentiation p itself could be made a time dependent variable, 
although this was not looked at here. The fractional derivative results are 
what might have been anticipated: however for consistency one should start 
with fractional derivatives at the first step and this leads to the subject of 
fractional derivative geometry. What fractional derivative geometry should 
look like is unclear, here a first attempt was made at guessing what curvature 
and line elements might look like in two dimensions. By chance this lead to 
a line element depending on the Gamma-function. 
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Figure 1: lowest line= 2.10^ exp(— r'^'^/O.l), middle line= 1/r, top line= 
5.10-5 exp(r-O-VO.l) 
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Figure 2: Red line is ^C^, blue points are y = 0.35(a; — 1) 
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